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Singularity of Some Random Continued Fractions 

BY Russell Lyons 



Abstract. We prove singularity of some distributions of random contin- 
ued fractions that correspond to iterated function systems with overlap and 
a parabolic point. These arose while studying the conductance of Galton- 
Watson trees. 



§1. Introduction. 

We are interested in the distributions of certain random continued fractions. They 
arose while studying the following question about Galton- Watson trees: Suppose that pk G 
[0, 1) for k > 1 with ^j^Pk = 1- Let T be the random genealogical tree resulting from the 
associated Galton- Watson process beginning with one individual, p, where each individual 
has k children with probability pk- Adjoin a parent A to p. Then the probability 7 that 
simple random walk starting at A will return to A is equal to the effective conductance of 
TU {A} from A to infinity when each edge has unit conductance (see, e.g., Doyle and Snell 
(1984) or Lyons and Peres (1998)). This quantity 7 enters, for example, in calculations of 
the Hausdorff dimension of harmonic measure on the boundary of T (see Lyons, Pemantle 
and Peres (1995)). Thus, it is of interest to calculate the distribution of 7. If denotes 
its cumulative distribution function, then satisfies the equation 



if s > 1 . 

Moreover, the functional equation (1.1) has exactly two solutions, and the Heaviside 
function l[o,oo) (see Lyons, Pemantle and Peres (1997)). Only a numerical computation of 
F-y is known. In all cases in which this calculation was carried out, it appeared that the 
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distribution was absolutely continuous with respect to Lebesgue measure with bounded 
density. Furthermore, unpublished work of the author and K. Zumbrun has shown that 
one can obtain a posteriori bounds on the modulus of continuity of Fry from such computer 
calculations. Therefore, it seems reasonable to conjecture, as in Lyons, Pemantle and Peres 
(1997), that the distribution of 7 is always absolutely continuous, at least when pfc = for 
all sufficiently large k. 

O. Haggstrom (personal communication, 1995) observed that one can write 

1 

7 = [1, 71, 1, 72, 1, 73, ...]:= 



1 + ^ 



1 

7i + 



1 + 

where 7^ are i.i.d. In the special case pi = P2 = 1/2, the random variables 7^ have the 
c.d.f. (l[o,oo) + F-y)/2. Note that [1, xi, 1, X2, . . . , 1, Xn] is monotonic increasing in each xf, 
this shows convergence of the infinite continued fraction. Thus, in attempting to prove 
the above conjecture on absolute continuity, we were led to consider the simpler random 
continued fractions 

[1, Xi, 1, X2, 1, X2, . . .] , 
where Xi are i.i.d. with the same distribution as 

X — i ^ with probability 1/2, 
I a with probability 1/2, 

where a e (0, 00) is a given real number. That is, instead of trying to solve a fixed point 
problem, we are attempting to understand the nature of the map taking a distribution 
of X to the distribution of [1, Xi, 1, X2, 1, X3, . . .] for a very simple X. Let /la be the 
distribution of this random continued fraction. 

Our main result is the following. Define etc to be the solution of Aq, = | log 2, where 
Aq. is defined in (2.6) below. 

Theorem 1.1. Ifa> 1/2, then is supported on a Cantor set of Hausdorff dimension 
< 1. If a < 1/2, then the support of is an interval. If a > ac, then is singular with 
respect to Lebesgue measure and is concentrated on a set of Hausdorff dimension < 1; we 
have ac G (0.2688, 0.2689). Whenever a > 0, the measure is continuous. 

The third sentence is the most novel one. 
Define 

_ , , X + a 
Ta{x) := 



1 + X + a 
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Thus, [l,xi,l , X2, • • • 5 1) Xn] — oTjjO- • -oTx^ (0). Note that Tq is not a strict contraction. 
The support of /la is contained in [0, Mq,] and includes the endpoints, where T^Ma = Mq, > 
0. If we solve this equation for Mq,, we get + aM^ — a = 0, so that 

—a + \J -\- ^OL 
Ma= - . 

In particular, M1/2 = 1/2 and T^O > TqMq, iff ck > 1/2. Since is supported on 
To[0, Ma] U Ta[0, Mq,], it follows by iteration that supp is a Cantor set for a > 1/2 and 
equals [0, Mq] for a < 1/2. Thus, Theorem 1.1 asserts that in some interval of a where 
there is overlap of To[0, Mq] and Tq[0, Mq], the measure //q is singular. This contrasts with 
the case of two linear maps, x ^ ax and x ^ a{l +x), where one has absolute continuity 
for almost all a in the overlap region (Solomyak 1995, Peres and Solomyak 1996). However, 
based on some numerical evidence, we believe that there is also an interval where //q is 
absolutely continuous: 

Conjecture 1.2. For all a sufficiently small, Ha is absolutely continuous with respect to 
Lebesgue measure. 

After seeing a preprint of this work, Simon, Solomyak, and Urbahksi (1998) made 
a great deal of progress on this conjecture. They proved that for Lebesgue-a.e. a e 
(0.215, CKc), the measure ji^ is absolutely continuous. In particular, this shows that the 
threshold etc in Theorem 1.1 is sharp. 

Other work on random continued fractions includes Bernadac (1993, 1995), Bhat- 
tacharya and Goswami (1998), Chamayou and Letac (1991), Chassaing, Letac, and Mora 
(1984), Kaijser (1983), Letac (1986), Letac and Seshadri (1995), Pitcher and Foster (1974), 
and Pincus (1983, 1985, 1994). In particular, Pincus (1983) asks about singularity of dis- 
tributions like Ha, when the images of two linear fractional maps overlap, which is our 
main contribution. Our analysis of the Hausdorff dimension is similar to that appear- 
ing on pp. 166-168 of Bougerol and Lacroix (1985), where they use Lyapunov exponent 
techniques to recover a result of Kinney and Pitcher (1965/1966). 

Finally, returning to our original problem about Galton- Watson trees, we mention 
another related problem: Show that Ylik>i absolutely continuous distribution, 

where is the size of the A;th generation of a Galton- Watson process. This arises as the 
resistance of the shorted Galton- Watson tree. In other words, if a graph is formed from a 
genealogical Galton- Watson tree by identifying all vertices in level k for each k separately, 
then ^X]fc>i -^fe^^) effective conductance from the original progenitor to infinity. 

Again, computation supports the conjecture that the distribution is absolutely continuous 
with a bounded density, at least when the offspring distribution is bounded. 
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In Section 2, we prove Theorem 1.1. In Section 3, we discuss the possibihty of or 
L°° densities for a < ac- 



§2. Proof of Theorem 1.1. 

By Theorem 6.5 of Urbahski (1996), the Hausdorff dimension of supp /la is less than 1 
when there is no overlap, i.e., when a > 1/2. (According to Remark 6.6 in that paper, we 
also have that the Hausdorff dimension of supp /la equals its upper Minkowski dimension.) 

Since fia is a stationary measure for the Markov chain on [0, Ma] that has transitions 
s I— > Txs, it follows that fi^ is continuous: By stationarity, for any s e [0, M^], we have 

In particular, if s is such that fiaii^}) i^ maximal, then ^^{{Tq^s}) — //q,({s}). Likewise, 
fia{{TQ~^s}) = fice{{s})- Since //q, is finite, it follows that //q,({s}) = 0. 
Define random variables An, Bn, Cn, -Dn by 

An Bn\ ^(1 Xi \ (1 Xn 
Cn Dn ■ V 1 l + Xj"'\l 1 + Xr. 



Since 



and 



1 X \ ( s\ _ ( s + Xt 
1 1 + X J \t J ~ \s + t + Xt 

s + Xt 1 



s + t + Xt ^ ^ _J_ 



we have that 



'^"'^^'^ =[l,Xi,l,X2,...,l,X„ + .], 



CnS + D 

n 

as is well known. We also write 

f an bn \ _ fl Xi \ fl Xn 

\Cn dn J ' \l 1 + Xi J \1 1+Xr 

for Xi = 0, a. 

Let y be a random variable with distribution fx^ independent of all Xj. Then 



fXa=mhXi,...,l,Xn + Y]) 
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for any n, where £(•) is the law, i.e., the distribution, of a random variable •. Now 

jC{[l,Xi,...,l,Xn + Y]) 

is a probability measure supported on the interval 

bn anMa + W 



[[1, Xi, . . . , 1, Xn], [1, Xi, . . . , 1, + Ma\\ = 



The length of this interval is Mo,/[dn{cnMo^ + dn)] since 



dn ' CnMa + dn 



^ndn b^Cfi — det 



Cn dn. 



Also, 

fXa[[[l,Xi, . . . ,l,Xn],[l,Xi, . . . ,l,Xn + Ma]]^ > 2"" , (2.1) 

with equality iff a > 1/2. By Billingsley's Theorem (Falconer (1997), p. 171), fia is 
concentrated on a set of Hausdorff dimension 



log 1///q ((x — e, a; + e)) 
^« -ess^sup liminf i^^yj^^ ' (2-2) 



which by (2.1) is at most 

lini log 2" _ log 2 ^2 3) 

log 

where Xa is the (top) Lyapunov exponent of the random matrix ^ j ^ ^ ^ ^ ; see (2.4) 

for the definition of Aq and Bougerol and Lacroix (1985), Cor. VI. 2. 3, for this property of 
Aq. [Simon, Solomyak and Urbahski (1998) have now proved that the quantity in (2.2) is 
in fact equal to that in (2.3) for Lebesgue-a.e. a G (etc, 1/2)-] In particular, is singular 
if Aq, > ^ log 2. Note that Aq, is strictly increasing in a since ^ j ^ +^^) stochastically 
increasing in a. Thus, if a > ac, then /Xq is singular and concentrated on a set of Hausdorff 
dimension less than 1. It remains to estimate ac- 

Parametrize projective one-space minus the horizontal direction by the vectors 

R X {1}. Leaving out the horizontal is of no consequence since ^ j ^ ^ ^ ^ maps it 
to another direction regardless of X. The random map ""^^ f J 1 ^ X ) ^ i^^^'^^^ ^ 
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random map (s, 1) i— > {Txs, 1) on R x {1} (thought of as P-^). The stationary measure for 
this Markov chain on R x {1} then becomes just /la itself. 
Now Aq, is the expected change in the fog norm on R^: 

A« = y"Efog[||(^J ^^^ys,l)\\/\\{s,l)\\]d^^M■ (2.4) 
If we use the norm 

||(a;,y)||e := ex + (1 - e)y 
for any choice of e e (0, 1), then we get 

f^'"(l, ea+(l-e)(l + a) , 1, e(a + a) + (1 - e)(l + g + a) ) ^ 
Jg^Q [2 es + (1 - e) 2 es + (1 - e) J 

(2.5) 

Since this holds for all e e (0, 1), it also holds for e = by the Bounded Convergence 
Theorem: 

Aa = -/ \og[{l + s){l + s + a)]diiais). (2.6) 

Define Ka to be the operator on c.d.f.'s on [0, Ma] corresponding to u t-^ (|To + |TQ,)zy. 
Thus, 

{KaF){s) = lF{T-h) + ^FiT-h) = lF{T-h) + ^FiT-'s - a) . 

Set Fo{s) := s/Ma (0 < s < Mq) and := KaFn- Then F„ converges to the c.d.f. 

F(q,) of fia because | [[1, xi, . . . ,1, Xn], [1, a^i, • • • , 1, + Ma]] \ < < 1/n. Since Ka is a 
monotone operator (i.e., if F < G, then KaF < KaG), we have that F(q,) < Fn+i < Fn 
whenever Fi < Fq. If a < 1/2, then Fi is linear from to Tq,0, from TqO to ToMq,, and 
from ToMa to Ma, with maximum slope in the middle portion. If a > 1/2, then Fi is 
linear from to ToMa, constant from TqMq to Tq,0, and linear from Tq,0 to Mq,. Thus, in 
both cases, it follows that Fi < Fq iff Fi(ToMa) < Fo{ToMa). Now 

= ^ = irk 

and 

1 a 

Fi{ToMa) = -[Fo{Ma) + Fo{Ma - a)] = 1 



2^^' , ^^,^.^a 2Ma 

Therefore, Fi < Fq iff a > 1/6. (We remark that M^/q = 1/3.) 

Now if F < G and h is an increasing function, then J h dF > J h dG. Since s 
log[(l + s)(l + s + a)] is increasing, this means that 



Aa > ^ y log[(l + s)(l + s + a)] dFn{s) 
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for all n, giving us a lower bound for each n; in particular, as computation shows, Aq, > 
I log 2 for a = 0.2689. Therefore, etc < 0.2689. 

We can get an upper bound on Aq, by a similar method. We have that (2.5) holds 
for e = 1 by the Bounded and Monotone Convergence Theorems (use monotonicity in a 
neighborhood of s = and boundedness elsewhere): 

1 /•^« / a\ 
Aa = 2 y log (^1 + -j (i^«(s) . 

(Note that for best numerical estimates, though without error bounds, one should choose 
a norm such that the expected log change is most nearly constant, rather than e = 0, 1 in 
(2.5). But the choices e = 0, 1 give us lower and upper bounds when cu > 1/6.) 

Now a F < G and h is a decreasing function, then J hdF < J hdG. Since s h- > 
log(l + a/s) is decreasing, we get the upper bound 

Aa < log (l + ^)(iF„(s) 

for all n; in particular, Aq, < | log 2 for a < 0.2688. This completes the estimate of ac- I 



§3. LP Densities. 

We next show that /la cannot be absolutely continuous with an density for a > 
V6/2 — 1 = 0.2247+. This result contrasts with the case of two equally likely linear 
maps (Solomyak 1995, Peres and Solomyak 1996), where all known instances of absolutely 
continuous measures have an L^ density. 

Proposition 3.1. If Ha is absolutely continuous with an L^ density, then a < \/6/2 — 1. 

Proof. Suppose that /la. has a density fa- Recall that Y has distribution //q. For xi, . . . ,Xn e 
{0, a}, the probability measure >C([1, xi, . . . ,1, Xn + Y]) has a density gn (depending on 
xi, . . . , Xn) supported on the interval 

hn anMg + hn 

_dn CnMa + dn _ 

Therefore, 

fa= Yl '^~''9n. 
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Considering only the diagonal terms and using the Cauchy-Schwarz inequality, we find 
that 



> 



E 2- 



•2n 



Xi .,. . . ^X-r, 



= E 2 



-2ndn{CnMo, + dn) 



DnjCnM^ + D^) 



Write 



Then 



J-Vin. • — I/O n 



L>2 = (o 1)(M„0M,) 








whence 



B[D^] = ( 1 ) E[M„ (g) M, 
By independence, we have 



/0\ 





Vi/ 



where 



R:=E 



1 X \ 

1 l+Xyl® 



1 X 
1 1 + X 



/I 0\ 
110 
10 10 

Vi 1 1 1/ 



/I 



+ 



2 * 



q: 

1 1 + q; a a + a 

1 a 1 + a a + a"^ 

\l 1 + a 1 + a (1 + a)^/ 



^1 a/2 a/2 a^/2 \ 

1 l + a/2 a/2 a/2 + aV2 

1 a/2 l + a/2 a/2 + a'^/2 

\l l + a/2 l + a/2 l + a + a'^/2j 
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The characteristic polynomial of R is 

t^t^ - {4: + 2a + a^/2)t^ + (6 + 4a + /2)t^ - (4 + 2a)t + 1 , 

which has its largest root > 2 iff a > \/6/2 — f . Clearly, the Perron- Frobenius eigenvector 
of R has all its coordinates strictly positive, so that 



^E[D2] ^ 2-^^ ( l)R' 






\l/ 



oo 



as n — > oo if a > -\/6/2 — 1. Hence //q, cannot have a density in for a > V6/2 — 1. I 

Similar methods show that if a > 3^/2 — 4 = 0.24"'", then cannot have a density in 
L^/^. In the other direction, we can extend the method to show the following proposition: 

Proposition 3.2. If is absolutely continuous with a density in for every p < oo, 
then a < (3-\/2 — 4)/2 = 0.1213"''. In particular, this is the case if fia has a density in L°°. 

Proof. Let p be a positive integer and 

Rp{a) := E[Mf '(^-')] , 

where Mf denotes the r-th tensor power ofMi:=(} Write r := 2(p- 1). 



1 1 + 

The method of proof of Proposition 3.1 shows that /Uq, cannot have a density in for 
ct > ctp, where ap is the value of a such that the largest eigenvalue 7p of Rp is equal to 
2(p-i)/2. Let the eigenvectors of Mi be ui and U2. Then the eigenvectors of 'M.f^ are 
■ ■ ■ ® Ui^, where all ij are either 1 or 2. Thus, the largest eigenvalue of M^^ is 5^^ 

1 /r 

where 6 is the largest eigenvalue of Mi. It follows that 7p tends to the largest eigenvalue 
of ^ J ^ ^ ^ ^ . It is then an easy matter to show that ap (3-\/2 — 4)/2. I 

Remark 3.3. Numerical evidence suggests that in fact may not have a bounded density 
until a is less than about 0.05 (very roughly). Furthermore, the density seems to gain 
increasing smoothness the smaller a becomes. 

Acknowledgement. I thank Grahame Bennett, Rick Kenyon, Yuval Peres, and Kevin 
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